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Abstract 


Payload  capabilities  were  calculated  for  an  expendable  upper 
stage  vehicle  compatible  with  the  Space  Shuttle  Vehicle.  Analysis 
was  performed  for  a four-stage  vehicle  that  was  modeled  with  impulsive 
thrust  and  transfer  trajectories  which  obey  restricted  two-body 
equations  of  motion. 

The  magnitude  of  the  maximum  payload  deployed  into  one  of  two 
specified  orbits  when  the  other  payload  is  known  is  solved  by 
breaking  the  four-impulse  transfer  into  two  dual-impulse  transfer 
trajectories.  The  maximum  payload  solution  for  one  transfer  depends 
upon  the  specified  payload  of  the  other  transfer.  Each  of  the  dual- 
impulse transfer  trajectories  is  determined  by  solving  a quartic 
equation  in  the  square  rbot  of  the  semi-latus  rectum  of  the  transfer 
orbit.  Maximum  payload  capability  was  dependent  upon  the  available 
impulse,  the  angle  between  terminal  orbit  planes,  the  difference  in 
the  radii  of  the  terminal  orbits,  the  plane  changes  at  departure  and 
arrival  points,  and  the  transfer  angle.  Transfer  solutions  were 
programmed  on  a CDC  6600  digital  computer. 

Computed  results  indicate  that  the  model  vehicle  is  capable  of 
many  non-coplanar  orbit-to-orbit  transfers  that  still  yield  practical 
payloads.  As  the  transfer  angle  deviates  from  the  neighborhood 
around  180°  and  the  other  geometrical  parameters  increase,  the 
payload  decreases. 


MAXIMUM  PAYLOAD,  FOUR- IMPULSE,  NON-COPLANAR, 
ORBITAL  TRANSFERS  FOR  AN  UPPER  STAGE  VEHICLE 
OF  THE  SPACE  TRANSPORTATION  SYSTEM 


I.  Introduction 


Background 

The  Space  Transportation  System  is  composed  of  a Space  Shuttle 
Vehicle,  a propulsive  upper  stage,  and  numerous  ground  systems 
necessary  for  support  functions.  The  Space  Shuttle  Vehicle, 
illustrated  in  Fig.  1,  consists  of  a payload-carrying  Orbiter  and 
three  external  propellant  tanks.  The  Space  Shuttle  Vehicle  will  take 
off  vertically  from  a launch  pad,  powered  by  the  three  mein  engines 
of  the  Orbiter  and  the  two  solid  rocket  boosters.  The  solid  rocket 
boosters  are  jettisoned  after  their  expiration.  Still  powered  by  its 
main  engines,  the  Orbiter  continues  ascent  into  earth  orbit.  The 
external  tank  is  jettisoned  just  prior  to  orbital  insertion  (Ref  2:419) 
Payloads  of  up  to  65,000  pounds  can  be  carried  in  the  cargo  bay,  which 
is  15  feet  in  diameter  by  60  feet  in  length.  The  onboard  maneuvering 
system  of  the  Orbiter  makes  it  possible  to  climb  to  an  orbital 
altitude  of  up  to  250  nautical  miles  (Ref  1:2.5. 1). 

The  propulsive  upper  stage  will  be  carried  into  low  earth  orbit 
in  the  paylaod  bay  of  the  Orbiter.  Once  on  station,  the  bay  doors 
will  open  and  the  upper  stage  vehicle  will  be  positioned  for  launch 
(see  Fig.  2).  DOD  has  selected  Burner  II,  an  expendable  interim  upper 
stage  (IUS)  vehicle,  to  perform  shuttle  missions  until  approximately 
1984,  at  which  time  NASA  plans  completion  of  a fully  reusable  upper 
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SHUTTLE  VEHICLE 

• DESIGN  LAUNCH 

WEIGHT 

AIMS  LBS. 

• FAYLOAO 

DUE  EAST 

GSK  LBS. 

DESIGN  |104‘| 

32S  LIS. 

0R8ITER 


PAYLOAD  UT 
IS  fl.  DM.  I (8  FT.  L0N6 

• 3 MIN  ENGINES 

470R  LBS  VAC  THRUST  EACH 

• BOMS  (HEINES 

EOQO  LIS.  VAC.  THRUST  EACH 

• 31  RCS  (HEINES 

S09  LBS  VAC  THRUST  EACH 

• B RCS  VERNIER  ENGINES 

25  LBS  VAC  THRUST  EACH 

• WEIGHT 

DRV  1SOK  LBS. 

DESIGN  LANDING  II7K  LBS. 


SOLID  ROCKET  BOOSTER 


• LENGTH  Ml  FT. 

• DIAMETER  MB  '* 

• WEIGHT 

LAUNCH  11S3K  LBS.  ‘A 
RECOVERY  I54R  LSS.  UtH 

• THRUST 

LAUNCH  2.SM  LBS.  EACH 


EXTERNAL  TANK 

• LENGTH 

153.1  FT. 

• OIAMETER 

330  1 IN. 

• WEIGHT 

LAUNCH 

IS31R  LBS. 

USABLE  PROP 

I550KLIS 

Fig.  1.  Space  S>uirrle  Vehicle 
(From  Ref  1 i: 35) 


Fig.  2.  Positioning  Upper  Stage  for  Launch 
(From  Ref  2:25) 


stage  vehicle. 

During  a typical  mission,  the  Space  Shuttle  Vehicle  delivers  the 
propulsive  upper  stage  and  satellite  into  orbit  around  the  earth.  The 
Space  Shuttle  Vehicle  releases  the  upper  stage  and  satellite,  as  a 
unit,  while  in  this  orbit.  The  propulsive  upper  stage  then  expends 
its  energy  to  maneuver  the  satellite  into  a different  orbit  around 
the  earth.  Eventually,  the  Space  Transportation  System  will  allow  for 
placement,  retrieval,  and  repair  of  orbiting  satellite  systems 
anywhere  around  the  earth  (Ref  7:14-17). 

Using  Burner  II  design  philosophy,  the  Boeing  Aerospace  Company, 
Space  Systems  Division,  has  developed  a family  of  solid  rocket  motor 
stages,  based  upon  a two-stage  model,  that  meets  the  requirements  for 
a low-cost,  expendable  interim  upper  stage.  Since  the  Burner  II  two- 
stage  vehicle  model  is  a short,  compact  unit,  two  will  fit  in  the 
Orbiter  payload  bay.  With  a four-stage  tandem  vehicle  configuration, 
it  is  possible  to  deploy  two  satellites  and  only  use  one  IUS  vehicle. 
As  a result,  the  shuttle  flights  needed  to  launch  spacecraft  that  use 
the  IUS  can  be  greatly  reduced.  For  this  reason,  this  thesis  analyzes 
a four-stage  tandem  vehicle  configuration.  Details  of  the  vehicle  are 
shown  in  Fig.  3.  The  weights  and  propulsion  summary  of  the  four-stage 
tandem  vehicle  configuration  under  analysis  appear  in  Table  I. 
Statement  of  the  Problem 

The  problem  is  to  calculate  the  maximum  transferrable  payload 
capability  of  a four-stage  solid  rocket  vehicle  that  is  capable  of 
deploying  two  payload  packages  into  two  different,  non-coplanar  orbits 
around  the  earth.  The  mass  of  one  payload  will  be  specified  so  that 
with  the  given  fixed-fuel  loading  the  unknown  payload  can  be  maximized 
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Payload  Payload 


Fig.  3.  Four-Stage  Vehicle  Model 


Table  I.  Model  Vehicle  Weights  and  Propulsion  Summary 


Total  Orbiter  Payload  Capacity 

65,000  LB 

Total  IUS  Weight 

54,700  LB 

Stage  1 

Stage  Weight 

20,400  LB 

Propellant  Weight 

18,000  LB 

Total  Inert  Weight 

2,400  LB 

Average  Thrust 

40,000  LB 

Average  Burn  Time 

135  SEC 

Stage  2 

Stage  Weight 

16,200  LB 

Propellant  Weight 

14,000  LB 

Total  Inert  Weight 

2,200  LB 

Average  Thrust 

40,000  LB 

Average  Burn  Time 

100  SEC 

Stage  3 

Stage  Weight 

12,000  LB 

Propellant  Weight 

10,000  LB 

Total  Inert  Weight 

2,000  LB 

Average  Thrust 

29,000  LB 

Average  Burn  Time 

100  SEC 

Stage  4 

Stage  Weight 

6,100  LB 

Propellant  Weight 

4,700  LB 

Total  Inert  Weight 

1,400  LB 

Average  Thrust 

14,700  LB 

Average  Burn  Time 

80  SEC 

with  respect  to  the  following  parameters: 

1.  The  total  impulse  available, 

2.  The  angle  between  the  planes  in  which  the  terminal  orbits  lie, 

3.  The  orbital  radius, 

4.  The  inclination  angles  required  to  transfer  from  the  initial 
orbit  into  the  transfer  orbit  and  from  the  transfer  orbit 
into  the  target  orbit,  and 

5.  The  transfer  angle  between  the  point  of  departure  and  the 
point  of  arrival  as  measured  between  the  radius  vectors. 

Only  trajectories  corresponding  to  the  shortest  transfer  arc  are 

considered  because  of  their  faster  transfer  times  and  usefulness  in 

orbital  rendezvous  missions.  The  satellite  or  other  payload  weight 

which  is  delivered  by  the  propulsive  upper  stage  can  vary  for  any 

particular  mission  as  long  as  the  combined  upper  stage  plus  payload 

does  not  exceed  65,000  pounds. 

The  purpose  of  this  study  is  to  investigate  the  maximum  payload 
capability  for  a proposed  tandem  configuration  of  a propulsive  upper 
stage.  Specifications  of  the  vehicle  being  analyzed  resemble  the 
Boeing  Company's  proposal  that  resulted  in  the  selection  of  Burner  II 
for  the  IUS  (Ref  10).  The  fundamental  assumption  was  that  the  vehicle 
employed  impulsive  thrust.  Other  assumptions  include  the  following: 

1.  The  vehicle  remains  in  orbits  about  a spherical  earth, 

2.  The  vehicle  is  free  from  perturbations  due  to  other  masses, 

3.  The  vehicle  is  free  from  drag, 

4.  The  vehicle  is  moving  under  the  influence  of  an  inverse 
square  gravitational  field, 

5.  No  gravitational  losses  occur,  and 
All  transfers  are  non-coplanar . 
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General  Approach 


(( 


The  approach  taken  to  solve  the  orblt-to-orbit  transfer  problem 
was  to  determine  the  maximum  payload  mass  for  which  a transfer 
trajectory  existed  between  points  In  two  terminal  orbits.  The  object 
was  then  to  define  the  boundaries  of  the  vehicle's  payload  capability 
by  extending  each  parameter  of  the  problem. 

The  four-impulse  trajectory  profile  is  a discontinuous  set  of 
two  dual-impulse  transfers  from  a known  point  in  the  initial  orbit  to 
a specified  point  in  the  target  orbit.  An  additional  constraint 
imposed  upon  the  problem  is  that  each  dual-impulse  trajectory  use  all 
of  the  total  impulse  available  from  the  two  solid  rocket  motors  required 
to  make  the  transfer. 

Each  dual-impulse  solution  is  in  the  form  of  a quartic  equation 
in  the  square  root  of  the  semi-latus  rectum  of  their  respective  transfer 
orbit.  The  coefficients  of  the  quartic  equation  are  functions  of  the 
geometry  of  the  transfer,  the  total  impulse  available,  and  the  velocity 
components  required  to  maintain  the  specified  terminal  orbits. 

Upon  factoring  the  quartic  equation,  candidate  values  for  the 
semi-latus  rectum  of  the  transfer  trajectory  result.  Maximum  payload 
is  defined  as  that  payload  mass  just  before  the  condition  that  no 
candidate  values  for  the  transfer  trajectory  exist.  No  candidate 
values  corresponds  to  the  point  where  the  total  system  mass  is 
increased  until  all  roots  of  the  quartic  equation  become  complex. 

Options  for  fixing  the  mass  of  the  first  payload  deployed  or 
fixing  the  mass  of  the  second  payload  deployed  are  incorporated  into 
the  computer  program.  The  maximum  payload  solution  then  depends 
upon  the  value  of  the  mass  specified  for  one  of  the  payloads.  The 
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major  resource  required  for  this  thesis  Is  the  CDC  6600  digital 
computer . 

Sequence  of  Presentation 

The  results  of  this  thesis  are  presented  as  follows:  In  Chapter 

II,  the  two-impulse  transfer  Is  analyzed.  A summary  of  the  development 
of  the  two-impulse  solution  is  presented.  Geometry  of  the  non-coplanar 
transfer  is  defined.  In  Chapter  III,  the  results  of  this  thesis  are 
presented.  Some  operational  limitations  for  the  model  analyzed  are 
defined.  The  thesis  is  concluded  with  a summary  of  conclusions  and 
recoamendations  for  further  study  in  Chapter  IV. 
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II . Calculation  of  Impulsive  Transfer* 


Impulsive  Thrust 

One  advantage  of  the  impulsive  thrust  approximation  is  that  it 
simplifies  what  might  otherwise  have  been  a complicated  analysis. 
However,  a disadvantage  is  that  it  is  not  extremely  accurate.  For  an 
impulsive  thrust,  the  magnitude  of  the  thrust  is  assumed  to  be 
infinite  and  applied  over  an  infinitesimally  small  time  period.  The 
impulsive  thrust  assumption,  therefore,  is  only  an  approximation  of 
the  more  realistic  program  of  finite  thrust.  However,  application  of 
the  calculus  of  variations  method  to  optimize  finite  thrust  trajectories 
generally  results  in  a nonlinear  two-point  boundary  value  problem. 

Closed  form  solutions  of  this  nonlinear  two-point  boundary  value 
problem  are  generally  not  available.  The  numerical  techniques  used 
to  solve  the  two-point  boundary  value  problem  often  require 
considerable  computer  resources,  time,  and  expense.  The  effect  of  an 
impulse  is  an  instantaneous  change  in  velocity  without  a corresponding 
change  in  position.  The  entire  transfer  trajectory  is  a conic  section 
which  obeys  the  restricted  two-body  equations  of  motion. 

In  Reference  13,  Rapp  has  shown  that  the  impulsive  solution  is  a 
good  first  approximation  to  the  finite  thrust  solution.  In 
Reference  8,  Robbins  presents  formulas  for  estimating  the  performance 
penalties  arising  from  the  use  of  impulses  instead  of  finite  thrust 
for  orbit  transfer  maneuvers.  Finally,  in  Reference  5,  Handelsman 
demonstrates  how  guessing  the  Lagrangian  multiplier  functions  A^(t) 
may  be  eliminated  by  use  of  the  Aj(0)  from  impulsive  trajectories 
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to  initiate  iteration  for  nonimpulsive  fixed-thrust,  fixed-exhaust 
velocity  propulsion.  The  impulsive  approximation  becomes  less  severe 
as  the  value  of  the  thrust  Increases  and  the  thrusting  periods 
decrease  with  respect  to  the  coasting  period  along  a thrust-coast- 
thrust  trajectory  profile. 

Maximum  Payload  Problem 

The  maximum  payload  is  constrained  by  the  total  amount  of  fuel 
available  for  the  transfer.  The  use  or  ejection  of  fuel  mass,  dm, 
results  in  an  increase  of  velocity,  dv.  If  the  velocity  and  initial 
mass  of  the  propulsive  upper  stage  at  some  time  t is  given  as  vQ  and 
mQ,  respectively,  and  the  exhaust  velocity  (measured  with  respect  to 
the  IUS  vehicle)  is  given  as  ve,  then  the  expression  for  the  change 
in  velocity  in.  an  inverse  square  gravitational  field  is 


AV  « V - v„ 


-Vglnfmo/m)  - g(r)dt 


If  the  exhaust  velocity  is  defined  as 


vg  =»  - thrust 
Am/  At 


and  gravity  losses  are  assumed  negligible,  then  the  change  in  velocity 
due  to  an  impulsive  thrust  becomes 


AV  » thrust  In (mo/m) 
Am/At 


If,  now,  u constraint  is  imposed  which  states  that  all  the  available 


fuel  be  used  in  t hi'  per forma nee  of  each  orbital  transfer,  then  the 


total  change  In  velocity  or  Impulse  becomes 

r - iT 


(4) 


where  AV^  ■ the  Individual  impulses  used  for  transfer  and 
IT  “ the  total  impulse  available. 

Analytical  expressions  for  the  calculation  of  the  change  in  velocity 
at  each  of  the  four  impulsive  thrusts  appear  in  Appendix  A. 

Coordinate  Frame 

In  order  to  generate  a solution  for  the  four-impulse  problem,  a 
suitable  coordinate  system  must  first  be  selected  to  describe  the 
initial  and  final  position  and  velocity  components  at  the  specified 
end  points.  In  this  non-coplanar  case,  a rotating  spherical  coordinate 
system,  illustrated  in  Fig.  4,  was  used.  The  symbol  U is  a unit  vector 

__  A 

in  the  radial  direction,  parallel  to  radius  vector,  R.  The  symbol  V 
is  a tangential  unit  vector  normal  to  U and  in  the  direction  of  the 

A 

orbital  motion.  The  symbol  W is  the  out  of  the  plane  component  of 
the  mutually  orthogonal  set  of  unit  vectors  that  is  defined  by  the 
cross  product  of  U and  V. 

Geometry  of  Impulsive  Transfer 

The  geometry  for  the  non-coplanar  two-impulse  transfer  is 
illustrated  in  Fig.  5.  The  vehicle  is  orbiting  the  earth  in  an 
initial  parking  orbit,  as  designated  by  orbit  1.  There  is  one 
instantaneous  thrust  from  the  initial  orbit  1 into  the  transfer  orbit, 
occurring  at  the  designated  initial  point.  The  change  in  velocity 
at  this  point  is  given  by. 
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A v s V — V 
n 1 vTj  V1 

There  is  a second  instantaneous  thrust  required  to  transfer  into  the 
final  orbit  2,  occurring  at  the  designated  target  point.  The  change 
in  velocity  at  this  point  is  given  by 


AV2  “ V2  - VT2 


(6) 


The  third  instantaneous  thrust  is  identical  to  the  first  thrust, 
except  it  transfers  the  vehicle  from  orbit  2 into  the  second  transfer 
orbit.  The  fourth  instantaneous  thrust  is  identical  to  the  second 
thrust,  except  it  transfers  the  vehicle  into  the  final  orbit  3.  The 
expressions  for  the  third  and  fourth  impulses  are 


and 


(7) 

(8) 


respectively.  The  four-impulse  trajectory  profile  is  a combination 
of  two  dual-impulse  trajectory  profiles.  For  example,  the  two- 
impulse  trajectory  profile  is  a direct  transfer  from  a known  point, 
designated  by  Rj  , to  a specified  point  in  the  target  orbit, 
designated  by  R£  • Another  direct  transfer  occurs  from  orbit  2 to 
orbit  3.  The  vehicle  must  deploy  the  first  payload  into  orbit  2,  so 
this  intermediate  orbit  is  also  the  orbit  in  which  any  necessary 
phasing  is  accomplished  between  maximum  payload  transfer  trajectories. 


IS 


This  would  seem  to  he  .1  rel.it  ivelv  simple  problem.  However,  when 

dealing  with  non-coplanar  orbits  and  when  not  transferring  from  an 
initial  position  along  the  intersection  of  the  orbit  planes,  each 
target  position  dictates  a different  total  amount  of  plane  change. 

This  plane  change  represents  a use  of  the  total  impulse  available. 
Geometry  of  Non-Coplanar  Transfer 

Other  geometry  for  the  non-coplanar,  orbit-to-orbit  transfer  is 
illustrated  in  Fig.  6.  The  arrows  designate  the  direction  of  orbital 
motion.  The  parameters  pictured  and  investigated  in  this  study  are 
defined  as  follows: 

1.  A,  the  angle  between  the  planes  in  which  the  terminal 
orbits  lie, 

2.  R,  the  radius  to  the  target  orbit, 

3.  ij,  the  inclination  angle  measured  from  the  initial 
orbit  to  the  transfer  orbit, 

4.  i2>  the  inclination  angle  measured  from  the  transfer 
orbit  to  the  final  orbit,  and 

5.  9,  the  transfer  angle  between  the  point  of  departure 
and  point  of  arrival. 

The  symbols  $ and  $ are  used  to  define  the  initial  position  and 
target  position,  respectively.  The  initial  and  target  positions  are 
referenced  to  the  line  of  intersection  of  the  orbit  planes  and  are 
positive  in  a direction  opposing  the  direction  of  the  orbital  motion. 
The  definitions  of  geometrical  parameters  in  this  problem  are 
important  in  the  development  which  follows  in  the  next  section. 

De vc 1 opment  of  a Two-Impulse  So  1 ut ion 

The  solution  of  the  four-impulse,  maximum  payload  problem  involves 
the  determination  of  two  transfer  conics  which  satisfy  the  position 
ind  velocity  boundary  conditions  of  the  terminal  orbits  and  maximize 
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Intersection  of 
Orbital  Planes 


Fig.  6.  Geometry  of  Non-Coplanar  Transfer 
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the  payload  for  each  duo  1 - impul so  transfer.  The  two-linpul so  transfer 

conic  solution  follows  the  same  development  as  that  used  by  Saxon 
(Ref  12).  Saxon  used  the  work  of  Gscobal  (Ref  4)  as  a guideline  in 
his  development.  The  solution  to  the  two- impulse  transfer  conic  is 
in  the  form  of  a fourth  order  polynomial  in  the  square  root  of  the 
semi-latus  rectum  of  the  transfer  orbit.  A summary  of  the 
development  of  the  polynomial  is  given  in  this  section.  A complete 
development  of  the  polynomial  may  be  found  in  Appendix  A of  Reference 
12. 

The  derivation  of  the  polynomial  is  based  on  relating  the 
transfer  velocity  vectors  at  the  initial  and  target  points  to  the 
semi-latus  rectum  of  the  transfer  orbit  while  meeting  end  conditions 
with  a given  change  in  velocity.  Using  the  coordinate  system  and 
geometry  defined  in  the  previous  sections,  the  position  vectors  for 
the  initial  and  target  points  are 

(9) 
(10) 
(ID 

Similarly,  the  velocity  vectors  for  the  initial  and  target  points  are 

(12) 

(13) 

(14) 

o 
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V1  - Vi + Vi 

V2  * Vx2U2  + VyzV2 

*3  * V3  + V3 


R1  “ R1U1 

A 

r2  * r2u2 


R3  " R3U3 


In  t lie  formal. ition  >>l  iliis  problem,  t lit  * cent  r i i i t v ami  (list. in  < from 


the  surface  of  the  earth  to  orbit  perigee  are  known.  Using  basic  two- 
body  orbital  mechanics,  the  vector  quantities  in  Eqs.  (12),  (13),  and 
(14)  are  calculated  directly. 

In  the  plane  of  the  transfer  conic,  the  equation  for  the  transfer 
velocity  vector  at  the  initial  point  is 


vTl  ■ *A, + VtAi  <is> 

where  U_  = radial  unit  vector 
*1 

V^.  » tangential  unit  vector  in  the  direction  of  orbit  motion 

Rj  * radius  to  initial  point 

fj  * angle  between  Rj  and  the  perigee  of  the  transfer  orbit 
The  general  conic  equation  in  polar  form  yields 


where  = transfer  arc  semi-latus  rectum 
e.j.  = transfer  arc  eccentricity 
0 = angle  between  Rj  and  R2 

R2  = radius  to  target  point 


(16) 

(17) 


Solving  for  from  Eqs.  (16)  and  (17);  and  after  some  other 

man  ipul.it  Ions  wc*  arrive  at  the  tangential  velocity  component 
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V 

R f - 1 

1 Ti  -rT 


(18) 


where  p *■  the  gravitational  constant. 
The  radial  velocity  component  is 


/ p/£.j  (1  - cos0)  / p£.^  (RjCosS  - Rj) 


(19) 


sin6 


R^R^sinB 


Substituting  the  radial  and  tangential  velocity  components  into 
Eq.  (15)  provides  the  transfer  velocity  vector,  V_  , at  the  initial 


L1 


point  as 


/Td  - CO80H 


T ✓ p (RjCosft  - RjHcp 


sin6 


RjI^sinO 


y~U  in 


(20) 


( 


In  a similar  manner,  the  transfer  velocity  vector 


V„  - R „U_  + R f V, 

T'2  2 T2  2 T2  T2 


takes  the  following  form: 
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• : * , ’ t ■ ■ *(  • ' 


_ /""u(cos0  - l)lT  /"p^R,  - R.cos6)lT 

VT  “ L ♦ uT 

2 sin0  Rjl^sinO  2 


4 

/T  iT 


R2  l2 


Vector  addition  on  the  geometry  as  illustrated  in  Fig.  5 gives 


AV1  " VTj  " V1 


AV2  " V2  * VT, 


By  definition,  since  we  are  restricted  to  the  plane  of  the  transfer 


conic. 


vi  - vui + <vvi 


v2  ■ (Vx2)U2  + <Vy2)V2 


The  constraint  that  the  total  impulse  available  be  entirely  used  can 


be  written 


2 1S  2 % 

(AVj)  + (av2)  - I 


Squaring  Eq . (26)  twice  and  simplifying  terms  yields 
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2 2 2 2 2 2 4 

(AVj  - AV2)  - 2IT(AV1  + AV2)  + IT  - 0 (27) 


2 2 

Expressions  for  AVj  and  AV2  are  obtained  by  taking  the  dot  products 


2 

AVX  - AVj  • AVX  (28) 

AV2  - AV2  • AV2  (29) 


2 2 

Substituting  for  AVj  and  AV2  from  Eqs.  (28)  and  (29)  into  Eq.  (27), 
multiplying  "by  , and  redefining  the  coefficients  of  powers  of 

yields 


4 3 2 

Ay  + By  + Cy  + Dy  + E - 0 (30) 


where  y ■ 

i 

The  coefficients  of  Eq.  (30)  are  defined  as  follows: 


2 2 

A - Bj  - 2ItB4 

2 

B - 26^3  - 2ItB5 

2 2 4 

C - 2BjB2  + B3  - 2I_j.Bg  + IT 

.2 

D - 2B2B3  - 2ItB6 
2 2 

E * B2  - 2ItB? 


Expressions  for  the  B values  appear  in  Appendix  B and  are  function* 
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of  tin?  following: 

1.  The  velocity  components  necessary  to  maintain  the  initial 
and  target  orbits, 

2.  The  plane  change  angles, 

3.  The  transfer  angle,  and 

4.  The  radius  of  both  the  initial  and  target  orbit. 

Solutions  to  the  transfer  from  the  initial  orbit  to  the 

intermediate  orbit  are  computed  using  one  polynomial.  Transfer  from 
the  intermediate  orbit  to  the  final  orbit  simply  requires  reimplemen- 
tation of  the  polynomial  expression  whose  derivation  is  summarized 
above.  The  fourth  order  polynomials  in  the  form  of  Eq.  (30)  are 
factored  using  a FORTRAN  extended  subroutine  called  DMULR.  DMULR 
finds  all  the  complex  roots  of  a polynomial  with  real  coefficients 
and  uses  double  precision  arguments.  Upon  factoring  the  quartlc 
equation,  four  values  for  the  semi-latus  rectum  of  candidate 
trajectories  resulted.  Normally  two  values  were  complex  while  the 
other  two  values  were  real.  If  too  much  impulse  was  available,  then 
all  four  roots  of  Eq.  (30)  could  be  real  and  all  four  transfer 
trajectories  were  possible  candidates.  If  insufficient  impulse  was 
available,  then  all  four  roots  became  complex  and  no  candidate  was 
assumed  to  exist  for  the  transfer  trajectory. 

Having  completed  the  transfer  solution  to  a given  orbit,  the 
total  vehicle  mass  was  then  increased  and  the  procedure  repeated. 

The  maximum  payload  mass  was  defined  as  the  mass  just  before  all  the 
roots  became  complex.  This  condition  corresponded  to  no  possible 
trajectory  for  the  increased  payload  mass  or,  equivalently,  no 
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possible  trajectory  for  the  available  impulse. 


With  a value  of  the  square  root  of  the  seml-latus  rectum,  , 
of  the  transfer  trajectory,  a transfer  velocity  vector  could  be 
evaluated.  Knowing  the  transfer  velocity  and  corresponding  position 
allowed  the  determination  of  all  the  classical  elements  of  the  transfer 
arc  by  means  of  the  equations  of  two-body  orbital  mechanics.  If  the 
coplanar  case  were  of  interest,  the  angles  ij  and  i 2 , which 

define  the  plane  changes  required  at  the  initial  and  target  points, 
are  set  equal  to  zero. 


III.  Results 


Normalized  System  of  Units 

To  simplify  the  computer  programming  and  the  input  and  output 

data,  a standard  set  of  geocentric  units  were  used.  The  unit  of 

distance,  1 DU,  is  defined  as  the  earth's  mean  equatorial  radius.  The 

unit  of  velocity,  1 DU/TU,  is  defined  by  the  velocity  of  an  object  in 

a circular  orbit  of  radius  1 DU.  For  this  set  of  units,  the  value  of 

3 2 

the  gravitational  constant,  u , is  1 DU  /TU  . The  conversion  from 
English  units  to  canonical  units  are  given  by 


-4 

1 

n.mi.  = 

2.903656x10 

DU 

-3 

1 

sec  ■ 

1.239444x10 

TU 

-5 

1 

ft/sec  “ 

3.855604x10 

DU/TU 

This  set  of  values  was  extracted  from  Appendix  A of  Reference  3. 
Comments  on  the  Computer  Program 

To  perform  the  four-impulse  transfer,  maximum  payload  transfers 
were  computed  in  two  steps.  The  first  transfer  was  from  the  initial 
to  the  intermediate  orbit.  The  second  transfer  was  from  the 
intermediate  orbit  to  the  final  orbit.  The  weight  of  the  vehicle 
before  departing  from  the  initial  orbit  included  four  solid  rocket 
motors  and  two  separate  payloads.  At  the  intermediate  orbit,  the 
first  payload  was  deployed.  The  first  and  second  solid  rocket  motors 
were  assumed  to  be  jettisoned  after  their  use.  The  weight  of  the 
vehicle  before  departing  from  the  intermediate  orbit,  therefore, 
included  the  remaining  two  solid  rocket  motors  plus  the  remaining 


25 


payload . 


Between  each  orbit  change  there  occurs  two  impulses  and  two  plane 
changes.  The  plane  change  at  the  initial  point  of  each  orbit  change 
was  through  an  angle  ij  and  at  the  final  point  was  through  an  angle 
i2  (see  Fig.  6).  Because  the  spherical  trigonometric  relationships 
held  true  only  for  orbits  of  the  same  altitude,  departure  and  arrival 
points  on  terminal  orbits  of  different  altitude  were  specified  by  the 
transfer  angle  only. 

Data  read  into  the  main  computer  program  included: 

1.  The  eccentricities  (e  = 0)  of  the  initial,  intermediate, 
and  final  orbits, 

2.  The  radial  distance  above  the  earth's  surface  of  each 
orbit, 

3.  The  average  thrust  for  each  of  the  four  solid  rocket 
motor  stages, 

4.  The  average  burn  time  for  each  of  the  four  solid  rocket 
motors,  and 

5.  The  stage  weight  (propellant,  structure,  and  motors) 
for  each  of  the  four  solid  rocket  motor  stages. 

In  addition,  one  payload  weight  and  the  orbit  in  which  it  was  to  be 

deployed  was  specified. 

In  subroutine  TRAJ,  the  eccentricity  and  serai-major  axis  of  the 
transfer  trajectory  were  computed.  Also,  the  total  impulse  available 
was  checked  for  each  possible  trajectory  to  insure  that  impulse 
constraints  were  satisfied. 

Payload  Definition 

In  this  thesis,  the  payload  weight  is  defined  to  be  the  useful 
payload  unloaded  at  the  destination  plus  any  weight  due  to  supporting 
structures,  such  as  a housing  container.  One  exception  is  made  to 
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tliis  attempt  U>  standardize  tin-  dot  in  i t ion  of  payload  throughout  the 

thesis.  This  exception  occurs  whenever  the  solution  to  the  maximum 
payload  is  sought  for  the  first  orbit  transfer  from  the  initial  orbit 
to  the  intermediate  orbit  and  the  payload  into  the  intermediate  orbit 
has  been  specified.  In  this  case,  the  payload  weight  includes  the 
third  and  fourth  stage  weights  plus  the  weight  of  the  payload  being 
carried  into  the  final  orbit.  Consequently,  this  defines  the  entire 
useful  payload  unloaded  at  the  intermediate  destination  and,  in  a way, 
is  consistent  with  the  original  definition. 

Maximum  Payload  Capabilities 

The  figures  and  tables  presented  in  this  section  are  all  obtained 
using  the  model  vehicle  weights  and  propulsion  specifications  that  are 
listed  in  Table  I.  Results  for  a vehicle  with  different  weights  and 
propulsion  specifications  are  presented  in  Appendix  C. 

Payload  values  presented  in  the  figures  and  tables  that  follow 
apply  to  the  capability  of  the  four  stage  IUS  after  it  has  been 
released  from  the  Orbiter.  If  the  IUS  and  the  two  payload  packages 
must  be  delivered  by  only  one  Space  Shuttle  Vehicle,  then  the  total 
payload  values  plus  IUS  weight  could  realistically  not  exceed  about 
62,500  pounds.  (A  remote  manipulator  arm  system  that  prepares  the 
IUS  for  launch  from  the  Orbiter  is  assumed  to  weigh  2,500  pounds.) 

Of  course,  larger  payloads  can  be  realized  for  some  transfers  if  the 
IUS  and  the  two  payload  packages  are  delivered  to  the  initial  orbit 
by  more  than  one  Space  Shuttle  Vehicle.  In  the  computation  of 
maximum  payloads,  five  parameters  were  considered  because  they  affect 
payload . 
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Annie  Between  Orbit  Planes.  Tlie  maximum  payload  trc>nd  that 

results  when  the  angle  between  the  planes  of  the  terminal  orbits 
increases  is  illustrated  in  Fig.  7.  The  maximum  deliverable  payload 
into  the  filial  orbit  decreased  steadily  as  the  angle  between  the  orbit 
planes  increased.  For  example,  a payload  of  4000  lbs  can  be  transferred 
to  a final  orbit  within  a 55°  plane  separation  of  the  intermediate 
orbit.  Parameter  values  for  a 10°  inclination  angle  from  intermediate 
orbit  to  final  orbit  appear  in  Table  II. 

It  is  of  interest  at  this  point  to  mention  the  results  that  have 
been  recorded  by  Porter  (Ref  6:640).  Porter  concludes  that  the  fixed 
two-impulse  system  can  deliver  the  design  payload  mass  to  a final 
orbit  within  a 30°  plane  separation  of  the  intermediate  orbit. 

Radius  to  Target  Orbit.  The  radius  of  the  intermediate  and  final 
orbits  were  varied.  For  a range  of  inclination  angles.  Fig.  8 
demonstrates  that  payloads,  for  example,  of  approximately  4000  lbs 
can  be  placed  into  synchronous  orbit.  But,  as  the  radius  to  the 
target  orbit  increases,  the  maximum  payload  decreases.  Also,  the 
maximum  payload  will  decrease  as  the  inclination  angles  increase. 

Table  III  shows  a decrease  in  the  maximum  payload  also  occurs  when 
there  is  a decrease  in  the  transfer  angle. 

Inclination  Angles.  The  determination  of  the  maximum  payload  is 
dependent  upon  the  total  plane  change  required  at  the  terminal  points 
of  the  initial  and  target  orbits.  The  relationship  between  the 
angles  of  these  plane  changes  and  payload  are  Illustrated  in  Figs.  9 
and  10.  Again,  payloads  of  4000  lbs  can  be  transferred  to  synchronous 
orbit  if  the  required  plane  change,  ij,  from  the  intermediate  orbit  to 
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Transfer  from  Intermediate  Orbit  to  Final  Orbit 


the  pi. me  of  the  transfer  trajectory  is  within  approximately  Jt>  . 

For  a range  of  values  of  ip  Fig.  10,  for  example,  shows  that  for  a 
plane  change,  i2»  within  approximately  40°,  the  model  vehicle  can 
deliver  a 4000  lb  payload. 

Transfer  Angle.  The  transfer  angle  between  departure  and  arrival 

points,  as  illustrated  in  Fig.  11,  is  another  factor  in  the  analysis 

of  maximum  payload  capabilities.  As  the  magnitude  of  the  transfer 

angle  departs  from  180°,  the  payload  decreases  rapidly.  The  maximum 

payload  also  decreases  when  the  vehicle  must  transfer  between  orbits 

requiring  increased  plane  changes.  In  addition.  Table  IV  demonstrates 

that,  again,  the  maximum  payload  is  found  to  decrease  as  the  radius  of 

the  target  orbit  increases.  A 5000  lb  payload,  for  Instance,  can  be 

delivered  to  synchronous  orbit  if  the  transfer  angle  is  within  ±10° 
o 

of  180  . But,  the  maximum  payload  decreases  to  3000  lbs  if  the 
transfer  angle  falls  within  ±40°. 


Transfer  from  Intermediate 


9.  Maximum  Payload  vs.  Inclination  Angle  into  Transfer  Orbit 
(Payload  »1  « 8000  lbs,  R,  - 2.8  DU,  R,  - 5.6  DU,  0 - 170°) 
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V.  Conclusions  and  Recommendations 


Summary  of  Conclusions 

Utilizing  an  expendable,  four-stage  vehicle,  this  thesis  has 
explored  the  maximum  payload  capabilities  of  one  possible 
configuration  of  a Burner  II  type  vehicle  to  perform  non-coplanar 
orbit-to-orbit  transfers.  The  maximum  payload  values  for  transport 
from  a low  earth  orbit  to  orbits  at  approximately  geosynchronous 
altitude  are  included.  This  range  includes  a majority  of  missions 
planned  for  utilization  of  the  Space  Transportation  System  during  the 
decade  of  the  1980s. 

An  impulsive  thrust  assumption  allowed  calculation  of  maximum 
payload  weights  to  be  obtained  for  a range  of  geometrical  parameters 
within  the  performance  capabilities  of  a four-stage  solid  propellant 
propulsive  upper  stage.  This  impulsive  thrust  assumption  allowed  the 
four-impulse,  maximum  payload  problem  to  be  considered  as  two  dual- 
impulse problems,  whose  solutions  were  in  the  form  of  two  quartic 
equations  in  the  square  root  of  the  semi-latus  rectum  of  the  transfer 
orbit.  Geometry  considerations  and  constraints  on  the  available 
impulse  permitted  investigation  of  maximum  payload  transfer  trajectories 
in  light  of  the  following  parameters: 

1.  The  angle  between  the  planes  in  which  the  terminal 
orbits  lie, 

2.  The  radii  of  the  terminal  orbits, 

3.  The  inclination  angle  from  the  initial  orbit  into  the 
transfer  orbit, 

4.  The  inclination  angle  from  the  transfer  orbit  into 


38 


i wit 


the  target  orbit,  and 


5.  The  transfer  angle  between  the  radius  vector  to  the  point 
of  departure  in  the  initial  orbit  and  the  radius  vector 
to  the  point  of  arrival  in  the  target  orbit. 

Results  of  the  parameter  study  reveal  an  additional  mission 
flexibility  beyond  that  of  a Hohmann-type  transfer  for  a solid  rocket 
motor  IUS  vehicle.  Many  non-Hohmann  orbit  transfers  are  within  the 
capabilities  of  vehicles  like  the  Boeing  Burner  II.  They  offer 
considerable  savings  in  mission  time.  Through  some  of  these  non- 
Hohmann  transfers  sufficient  payload  mass  can  be  delivered  to  make 
this  vehicle  a worthy  choice  for  earth  orbital  missions  using  the 
Space  Transportation  System. 

Recommendations  for  Further  Study 

The  results  obtained  in  this  study  can  be  extended  in  the 
following  manner: 

1.  Use  the  more  realistic  theory  for  finite  thrust  to  set  up  a 
nonlinear  two-point  boundary  value  problem.  By  numerical  methods, 
the  solution  will  result  in  more  realistic  values  of  the  maximum 
transferrable  payload. 

2.  Investigate  both  impulsive  and  finite  thrust  solutions  for 
orbits  that  are  not  circular  but  exhibit  some  eccentricity. 
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Appendix  A 


Expressions  for  the  Calculation  of  Impulses 


If  the  total  vehicle  system  mass,  M,  is  defined  as 


M = Ml  + M2  + PAYLOAD  //I  + M3  + M4  + PAYLOAD  #2 

where 

Ml  = (Ml)propellant  + (Ml) inert 
M2  = (M2)propellant  + (M2) inert 
M3  = (M3) propellant  + (M3) inert 
M4  = (M4)propellant  + (M4) inert 


then,  the  impulses  for  each  stage  are  given  by  the  following: 


App, ■ mlix  B 


Expressions  for  Coefficients  of  the  Quartic  Equation 


where  y * fc-j. 


Ay4  + By3  + Cy2  + Oy  + E - 0 
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h 
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Append  i x C 

Plots  of  Computer  Results  for  a Second  Vehicle 


Table  V.  Another  Model  Vehicle's  Weights  and  Propulsion  Summary 


Total  Orbiter  Payload  Capacity 

65,000 

LB 

Total  IUS  Weight 

51,000 

LB 

Stage  1 

Stage  Weight 

24,000 

LB 

Propellant  Weight 

20,000 

LB 

Total  Inert  Weight 

4,000 

LB 

Average  Thrust 

42,000 

LB 

Average  Burn  Time 

140 

SEC 

Stage  2 

Stage  Weight 

12,000 

LB 

Propellant  Weight 

9,500 

LB 

Total  Inert  Weight 

2,500 

LB 

Average  Thrust 

14,000 

LB 

Average  Burn  Time 

100 

SEC 

Stage  3 

Stage  Weight 

12,000 

LB 

Propellant  Weight 

9,500 

LB 

Total  Inert  Weight 

2,500 

LB 

Average  Thrust 

14,000 

LB 

Average  Burn  Time 

100 

SEC 

Stage  4 

Stage  Weight 

3,000 

LB 

Propellant  Weight 

2,300 

LB 

Total  Inert  Weight 

700 

LB 

Average  Thrust 

6,000 

LB 

Average  Burn  Time 

60 

SEC 

Transfer  from  Intermediate 


o 
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Fig.  12.  Maximum  Payload  vs.  Angle  Between  Orbit  Planes 
(Payload  #1  - 1000  lbs,  R,  - R~  - 4.0  DU,  <|>  = 170°) 
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from  Intermediate 


Fig.  14.  Maximum  Payload  vs.  Inclination  Angle  into  Transfer  Orbit 
(Payload  #1  - 8000  lbs,  R,  - 2.8  DU,  R,  - 5.6  DU,  0 - 170°) 


Transfer  from  Intermediate 


Fig.  15.  Maximum  Payload  vs.  Inclination  Angle  into  Final  Orbit 
(Payload  #1  - 4000  lbs,  R - 2.8  DU,  R - 5.6  DU,  6 - 170°) 


Transfer  from  Intermediate 


90.00  135.00  • 180.00  225.00  270.00  315.00  360.00 

TRANSFER  ANGLE  (DEG) 


Table  VIII.  Maximum  Payload  vs.  Transfer  Angle 
Transfer  from  Intermediate  Orbit  to  Final  Orbit 
(Payload  #1  * 1000  lbs,  R7  = 1.5  DU,  i.  =i7  * 2 


of  payload  values  exist  for  9 >180' 


Append l x D 


Listing  of  Computer  Program 

In  the  pages  of  this  appendix  appear  the  computer  program  that 
was  used  to  obtain  data  for  this  thesis.  In  order  that  reading  these 
pages  might  be  alittle  easier  the  following  synonyms  are  defined: 


Ml 

M2 

M3 

M4 

M 

LOAD 

PAY 

TH1 

TH2 

TH3 

TH4 

BT1 

BT2 

BT3 

BT4 

PDAES1 

PDAES2 

PDAES3 

EX1 

EX2 

EX3 

EXT 

XINC1 

XINC2 

THETA 

SMAXIS 

TIMP 

X IMP /SIMP 

M1P 

Mil 

M2P 

M2 1 

M3P 

M3I 

M4P 

M4I 


stage  weight  of  SRM  It  1 (lbs) 
stage  weight  of  SRM  It 2 (lbs) 
stage  weight  of  SRM  It 3 (lbs) 
stage  weight  of  SRM  It 4 (lbs) 
total  upper  stage  vehicle  weight 

1 means  payload  to  be  unloaded  into  intermediate  orbit 

2 means  payload  to  be  unloaded  into  final  orbit 
designated  payload  mass 

average  thrust  of  SRM  #1  (lbs) 

average  thrust  of  SRM  It 2 (lbs) 

average  thrust  of  SRM  It 3 (lbs) 

average  thrust  of  SRM  #4  (lbs) 

average  burn  time  of  SRM  It  1 (sec) 

average  burn  time  of  SRM  #2  (sec) 

average  burn  time  of  SRM  It 3 (sec) 

average  burn  time  of  SRM  It 4 (sec) 

parking  orbit  distance  above  earth's  surface  (DU) 

intermediate  orbit  distance  above  earth's  surface  (DU) 

final  orbit  distance  above  earth's  surface  (DU) 

eccentricity  of  parking  orbit 

eccentricity  of  intermediate  orbit 

eccentricity  of  final  orbit 

eccentricity  of  transfer  orbit 

inclination  angle  between  initial  orbit  and  transfer  arc 

inclination  angle  between  transfer  arc  and  target  orbit 

angle  between  position  vectors 

semi-major  axis 

total  impulse  used 

total  impulse  available 

propellant  weight  of  SRM  //I  (lbs) 

inert  weight  of  SRM  It  1 (lbs) 

propellant  weight  of  SRM  It 2 (lbs) 

inert  weight  of  SRM  It 2 (lbs) 

propellant  weight  of  SRM  //3  (lbs) 

inert  weight  of  SRM  II 3 (lbs) 

propellant  weight  of  SRM  //4  (lbs) 

inert  weight  of  SRM  It 4 (lbs) 
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PROGRAM  MAIN  (INPUT*  OUTPUT  * TA  PE  6=0 JT PUT  »PLOT) 

C OMM ON 7P OUT /W 1, W2, W3, W4, W5, AL PM A1,ALPHA2,ALPHA3, ALPHA 4, ALPHAS, 
iAL»-»  A6,  ALPMA7,  ALPHAS 
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END 
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* AL°HA5 

OV2SO  = (ALPHA1*P)  ♦ (ALPHA6*SQRP>  ♦ (AL«»HA7/SQRP ) ♦ (ALPHA4/P) 
♦ A. PH AS 


0V1  * SORT (OV1SQ) 
OV2  = SORT ( OV2SD 
TIM3  s DV1  ♦ 0V2 
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PRD3  RAM  X PLOT  (INPUT,  OUTPUT,  HP*  5=1 MPUT, PLOT ) 
DIMENSION  X (22) ,Y (4,22)  ,YPLOT (22) 

INTEGER  ORBIT 

RENO (5 , * ) ORBIT ,LOAD,N»M 
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CM..  SYMBOL  (2.2, 5.4,  0.133,  19HSED0N0  ORBIT  CHANGE , 0.,  19) 
CALL  SYMBOL (2. 2, 5.2, 0. 150, 17H1ST  PAYLOAD  KNOWN, 0. ,17) 

CAL-  PLOT  (3.,  0., -3) 

CALL  PLOTE 
END 


sua*ojrrNF  am5l?( xiMCi,PHr , a, xr mc?, theta > 

°I  = ?.141?9265 
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The  maximum  payload  solution  for  one  transfer  depends  upon  the  specified 
payload  of  the  other  transfer.  Each  of  the  dual-impulse  transfer 
trajectories  is  determined  by  solving  a quartic  equation  in  the  square  root 
of  the  semi-latus  rectum  of  the  transfer  orbit.  Maximum  payload  capability 
was  dependent  upon  the  available  impulse,  the  angle  between  orbit  planes, 
the  difference  in  the  radii  of  the  terminal  orbits,  the  plane  changes  at 
departure  and  arrival  points,  and  the  transfer  angle.  Transfer  solutions 
were  programmed  on  a CDC  6600  digital  computer. 

Computed  results  indicate  that  the  model  vehicle  is  capable  of  many 
non-coplanar  orbit-to-orbit  transfers  that  still  yield  practical  payloads . 
As  the  transfer  angle  deviates  from  the  neighborhood  around  lSO^^n^  the 
other  geometrical  parameters  increase,  the  payload  decreases.  ' At  r> 
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